Analysis of the Airport Network of India as a complex weighted network 
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Transportation infrastructure of a country is one of the most important indicators of its economic 
growth. Here we study the Airport Network of India (ANI), which represents India's domestic civil 
aviation infrastructure, as a complex network. We find that ANI, a network of domestic airports 
connected by air links, is a small-world network characterized by a truncated power-law degree 
distribution, and has a signature of hierarchy. We investigate ANI as a weighted network to explore 
its various properties and compare them with their topological counterparts. The traffic in ANI, as 
in the World-wide Airport Network (WAN), is found to be accumulated on interconnected groups 
of airports and is concentrated between large airports. In contrast to WAN, ANI is found to be 
having disassortative mixing which is offset by the traffic dynamics. The analysis indicates toward 
possible mechanism of formation of a national transportation network, which is different from that 
on a global scale. 
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I. INTRODUCTION 

Transportation infrastructures are of crucial impor- 
tance to the development of a country [l[ and are im- 
portant indicators of its economic growth. They form 
the backbone of tourism industry, support movement of 
goods and people across the country, thereby driving the 
national economy. Roadways, railways, and airways are 
the major means of transport in India, although con- 
tribution of airways is small compared to that of the 
other two. The civil aviation sector in India has been 
developing steadily and is expected to grow at an ac- 
celerated pace with the shift in the policy of the gov- 
ernment, and the addition of several low-cost private air 
service providers Q . These private players offer compet- 
itive pricing, provide more region-oriented service, and 
thereby help to increase the air traffic. Understanding of 
these transportation systems is important for reasons of 
policy, administration and efficiency. 

During the past few years, complex network analy- 
sis has been used to study transportation systems (rail- 
ways i, airlines @, |, i 0, i), which are man-made 
infrastructures, from different aspects, apart from many 
other systems of diverse origins [3, [l(| • The World-wide 
Airport Network (WAN) has been studied from topolog- 
ical as well as traffic dynamics perspective. WAN has 
been studied [|[ for its degree distribution and a model 
with constraints (such as cost of adding links to the nodes 
or the limited capacity of nodes) was proposed to account 
for the truncation for high-degrees in its scale-free cumu- 
lative degree distribution. It has been observed [j| that 
in WAN, most connected nodes are not necessarily the 
most "central" - nodes through which most of the short- 
est paths go. A model, incorporating the geo-political 
constraints, has been proposed Q to explain this appar- 



ently surprising result. Beyond the topological features, 
WAN was studied Q in more detail by considering the 
traffic - strength of interactions between nodes - dynam- 
ics on it. A model was proposed [ll| for the evolution of 
weighted evolving networks in an effort to understand the 
statistical properties of real-world systems. Airport Net- 
work of China (ANC) 0, a network much smaller than 
WAN, was also analyzed for its topology and traffic dy- 
namics. Its topology was found to be having small-world 
network features and a two-regime power-law degree dis- 
tribution. 

We investigate India's domestic airport network which 
comprises air services of all major civil air service 
providers. We study the network for its topological fea- 
tures and for its traffic dynamics, by considering the in- 
tensity of interactions. First, we study the network as an 
"unweighted network" to investigate its topological prop- 
erties. This network is formed by considering whether or 
not a pair of airports are connected by an airline. As in 
many other complex networks, the details of flow of in- 
formation is a crucial factor for airport networks. Hence 
we study "weighted network" , in which we consider the 
number of flights plying between any two airports per 
week. "Unweighted network" represents the architecture 
(topology) of airports' connectivity, whereas "weighted 
network" represents traffic dynamics in the network. 

Our analysis shows that while ANI is similar to the 
WAN in some aspects, it has differences in some features 
as reflected in its network parameters. We find that ANI 
has small-world [12| network features and has a scale- 
free [HI degree distribution. The traffic in ANI is found 
to be accumulated on interconnected groups of airports 
and is concentrated on trunk routes between large air- 
ports. It is found to be having a disassortative topology, 
in contrast to WAN. 



II. AIRPORT NETWORK OF INDIA (ANI) 
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The Airport Network of India (ANI) comprises domes- 
tic airports of India and airlines connecting them. The 
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air services provided in ANI include all major domestic 
air service providers in India, including some interna- 
tional airlines which ply domestic routes. The data was 
obtained from a timetable (EXCEL 's Timetable for air 
services within India [23|, as on 12th January 2004) of 
air services, which includes domestic services provided by 
following air-travel service providers: Indian Airlines, Al- 
liance Air, Jet Airways, Air Sahara, Air Deccan, Jagson, 
Druk Air, Air India, Bangladesh Biman, Royal Nepal and 
Srilankan Airlines. 



A. Unweighted ANI 

The "unweighted ANI" is a directed network with 
(iV=) 79 nodes (airports) and 442 directed links (flights 
going from one airport to another). It is represented 
by a binary adjacency matrix, A(N x N), whose ele- 
ments dij take value 1 if there is a flight by any service 
provider from airport i to airport j on any day of the 
week and otherwise. The asymmetric adjacency matrix 
(A) was used to find properties (in-degrees, out-degrees, 
and shortest paths) which are sensitive to the direction. 
Certain properties (degrees, degree-correlations, cluster- 
ing coefficients) were calculated after symmetrizing the 
adjacency matrix, which is justified, as out of total 228 
flight routes 221 are bi-directional. We needed to add 
only 14 fictitious flights so as to symmetrize A. The to- 
tal number of edges in the symmetrized adjacency matrix 
is, M = 228. 



B. Weighted ANI 

To include the information about the amount of traffic 
flowing on the network, the "weighted ANI" is defined by 
considering the strengths of the links in terms of number 
of flights per week. It is represented by a weight matrix, 
W, where each element Wij stands for the the total num- 
ber of flights per week from airport i to airport j. Since 
majority of the nodes are symmetric in ANI, in terms of 
traffic flowing into and out of it (wij = Wji), we sym- 
metrize the weight matrix W which is used for all the 
weighted analyses. 



for the number of flights terminating-into and number 
of flights originating-from that airport, respectively. We 
observe that for a very large number of nodes in ANI 
kj n = k° ut . The average degree of symmetrized ANI was 
found to be (K) = 2M/N = 5.77. The average shortest 
path length (L) for a directed network with N nodes is 
defined as, 



L = 
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N(N — 1) 
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where = shortest path length from node i to j. Clus- 
tering coefficient (Ci) of a node is defined as the ratio of 
number of links shared by its neighboring nodes to the 
maximum number of possible links among them. In other 
words, Ci is the probability that two nodes are linked to 
each other given that they are both connected to i. The 
average clustering coefficient is defined as, 
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where Ci = clustering coefficient of node i. 

Small-world [H, G3] networks are characterized by a 
very small average shortest path length (L) and a high 
average clustering coefficient (C). Shortest path length 
from node i to j, Lij, is the number of flights needed 
to be taken to go from i to j by the shortest route. We 
found the average shortest path length of ANI to be (L-) 
2.2593, which is of the order of that of a random network 
{Lrand ~ hxN/ In (K) — 2.493) of same size and average 
degree. The average clustering coefficient of ANI was 
found to be (C=) 0.6574, which is an order of magni- 
tude higher than that of the comparable random network 
(Crand ~ (K) / N = 0.0731). These two properties indi- 
cate that ANI is a small- world network. WAN [1,0] and 
ANC [6] have also been found to be small-world networks. 
In particular, ANC, which is of comparable size of that 
of ANI, has average shortest path length and clustering 
coefficient, 2.067 and 0.733, respectively. 



A. Shortest Paths Analysis 



III. TOPOLOGICAL ANALYSIS OF 
UNWEIGHTED ANI 

Degree of a node is the number of nodes it is directly 
connected to. Degree of a node i is defined as, 

N 

i=i 

In a directed network, in-degree (out-degree) of a node 
is the number of in-coming (out-going) links. In ANI, 
in-degree and out-degree (k° ut ) of an airport stand 



ANI (a directed connected network with 79 nodes) has 
6162 distinct paths - node-to-node flight routes. Shortest 
path analysis is intended to give us an idea about the ease 
of travel in the network. Fig. [Tj shows the distribution of 
the shortest paths in ANI. 

Table Q] summarizes the results of shortest paths anal- 
ysis. Apart from other statistics, we also show number 
of flights needed to be changed, which is an indicator of 
the convenience of travel in the network; the lesser the 
better. Around 99% of paths are reachable by changing 
a maximum of 2 flights. Around 68% of paths are reach- 
able by changing maximum 1 flight. Whereas about 7% 
paths are connected by direct flights. We also deduce 
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FIG. 1: Shortest path distribution in ANI 



TABLE I: Break-up of the number of flight routes having 
a certain shortest path, their percentage and corresponding 
number of flights needed to be changed by the shortest route. 
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Percentage 


No. of Flights 


Path 


Paths 


of flight 
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routes 


changed 
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442 


7.1720 
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3741 


60.7100 


1 


3 


1918 


31.1260 


2 


4 


61 


0.9899 
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that the diameter (defined as the longest of all shortest 
paths) of ANI is 4, which means that one needs to change 
at most 3 flights to reach from any airport to any other 
airport in ANI. 

From the passengers' point of view, best is as few 
change of flights as possible, to reach from any airport 
to any other airport in the network. Though this is an 
ideal situation for the passengers, it is not economically 
viable for the air service providers. We infer that the 
airport network has evolved, with the efforts from the 
service providers, to cater for the convenience of passen- 
gers, thus acquiring the small-world topology. 

B. Degree Distribution 

Degree, as defined in Eq. [H is one of the measures of 
centrality of a node in the network. Degree symbolizes 
the importance of a node in the network - the larger 
the degree, the more important it is. The distribution 
of degrees in a network is an important feature which 
reflects the topology of the network. It may shed light 
on the process by which the network has come into ex- 
istence. The networks in which the links between two 
nodes are assigned randomly have a Poisson degree dis- 
tribution fl5j with most of the nodes having a typical 
degree. A large number of networks ranging from In- 




k (Degree) 

FIG. 2: Cumulative degree distribution in ANI. It can be 
approximated by the power-law fit, P(k) ~ k 1cum . 



ternet to protein-protein interaction network in yeast, 
have been found Q to be having scale-free degree dis- 
tributions. The scale- free distribution, characterized by 
a power law - P(k) ~ fc -7 , with a scaling exponent 7 
- can be explained with the help of a growing network 
model with preferential attachment of the nodes which 
are being added to the network jl3| . 

Since ANI is a small network, we analyze the cumula- 
tive degree distribution, P(> k), whose scaling exponent 
leum is related to that of P(k) by 7 = 7 cum + 1. We find 
that the cumulative degree distribution of ANI follows a 
power law as seen in Fig.[5]for a wide range of k, although 
it deviates from it for large degrees. The deviation for 
large degrees can be attributed to the cost of adding links 
to nodes or to the limited capacity of nodes. With the 
help of numerical simulations, it was shown Q that the 
cost of adding links leads to a cutoff of the power-law 
regime for large k in the cumulative degree distribution, 
as is the case with ANI. We find the scaling exponent of 
ANI to be, 7 = 2.2 ±0.1. 



IV. CENTRALITY IN WEIGHTED ANI 

We analyzed the weighted ANI, by considering the 
flow of information (traffic) on the topology of the net- 
work. The statistical analysis of weights between pairs of 
nodes indicates the presence of right-skewed distribution 
as shown in Fig.[3J This shows a high level of heterogene- 
ity in ANI, as also found in the case of Airport Network 
of China Q and WAN Q. 

It has been observed that the individual weights do 
not provide a general picture of the network's complex- 
ity [Ty] . We use a different measure of importance of a 
node which considers the flow of information in a net- 
work, following Barrat et al. The weighted counter- 
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FIG. 3: Cumulative weights distribution in ANI. 



FIG. 4: Average strength s(k) as a function of degree (k) of 
nodes. 



part of degree, strength (sj), is defined as, 

N 
3=1 

Strength of an airport represents total traffic handled by 
it per week. Fig. |4] shows the correlation between degree 
(fc) and the average strength of vertices with degree A;, 
s(k). We find that s(k) increases with k as, s(k) ~ k@ . 
If the strength and degree of a node were uncorrelated, 
then Si = (w)ki, where (w) = (2M) -1 J2i j a ij w iji which 
will yield (3=1. This is an uninteresting case, as in 
this situation weights do not provide any better infor- 
mation than degrees. For ANI, we find that the expo- 
nent to be (3 an i — 1-43 ± 0.06. This implies that the 
strengths of nodes are strongly correlated to their de- 
grees in ANI. Larger is an airport, the more traffic it 
handles, in contrast to what would be expected if the 
topology and the traffic were uncorrelated. This feature 
of ANI is similar to that of WAN, where it was found 
that (3wan — 1-5 ± 0.1. 



V. STRUCTURAL ORGANIZATION IN 
WEIGHTED ANI 

Various properties of network are decided, among 
many other things, by the way the network has evolved, 
the inherent interdependencies of the nodes, and the ar- 
chitectural constraints. Network parameters are defined 
so as to be able to capture such features of networks. We 
study few such parameters related to the local cohesive- 
ness and connection tendency of nodes, with and without 
consideration of weights to gain better insight of traffic 
dynamics on the architecture of ANI. 



A. Clustering Coefficient 

Clustering coefficient as defined in section [TTTl captures 
the local cohesiveness of a node. Average clustering coef- 
ficient (C) measures the global density of interconnected 
nodes in the network. Additional information at the in- 
termediate level can be obtained by defining C(k), aver- 
age clustering coefficient of nodes with degree k. These 
topological parameters overlook the flow of information 
on the network and hence may not present correct infor- 
mation about the network dynamics. The clustering co- 
efficient could be redefined [7| , to incorporate the weights 
of the edges, as, 

w I sr^Wij+Wih 

Ci = ^—1) ^ (5) 

This parameter measures local cohesiveness, by taking 
into account the interaction intensity found on the lo- 
cal triplets. The unweighted and weighted clustering co- 
efficients could be compared (the normalization ensures 
that < cf < 1) to assess the tendency of accumulation 
of traffic on interconnected triplets. Weighted cluster- 
ing coefficients averaged over all nodes (C w ) and over all 
nodes with degree k (C w (k)) are defined analogous to 
their topological counterparts. 

It has been shown [l7[ that hierarchical networks are 
expected to have a non-trivial, power-law decay of C(k) 
as a function of k, which means low degree nodes be- 
long to interconnected communities. Many real networks 
have been found [13, [H, [HI to be having such nontriv- 
ial decay. We find that in ANI, C(k) shows a power-law 
decay with C(k) ~ k . This points toward an inherent 
hierarchy in the architecture of ANI. This is an inter- 
esting result as it has been found [l7| that some other 
infrastructure networks (Internet at the router level and 
power grid network of Western United States) which are 
constrained by "geographical organization" do not show 
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FIG. 5: Average unweighted (C(fc)) and weighted (C™(fc)) 
clustering coefficients of nodes with degree fc. 
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FIG. 6: Average unweighted (k nn (k)) and weighted 
degree of nearest neighbors of nodes with degree k. 



such hierarchy. This could be because, as opposed to 
the above-mentioned networks, where the cost of having 
a link grows in proportion to the distance between two 
nodes, the geographical constraint is not so strong in the 
case of airport networks. It is relatively easier to add an 
air connection to the existing network, than to add a data 
cable or an electrical cable for comparable geographical 
distances. 

As seen in Fig. [5J after remaining almost constant for 
k < 8, C(k) falls rapidly thereafter. This indicates that 
large airports provide air connectivity to far-off airports, 
which themselves do not tend to be connected. In ANI, 
regional and national hubs provide air connectivity to 
airports in their domains, with the latter providing it to 
a much larger number, which reduces the C(k) of these 
hubs. We find that C w /C S 1.075, which means that 
the traffic is accumulated on interconnected groups of 
nodes which form high traffic corridors, also known as 
trunk routes. We find that C w (k) is restricted in its range 
across k and is consistently higher than the corresponding 
C(k). The higher the k, more pronounced is the differ- 
ence between C w (k) and C(k). This implies that with 
increasing degree, airports have progressive tendency to 
form interconnected g roup s with high-traffic links. This 
rich-club phenomena [20(, in which high degree nodes 
tend to form cliques with nodes with equal or higher de- 
gree, is also observed in WAN Q- 



B. Degree Correlations 

Another parameter which is used to investigate the 
networks' architecture is the average de gree of nearest 
neighbors, k nn (k), for nodes of degree k [2lj |. It is con- 
venient to calculate k nn (k), which is related to the cor- 
relations between the degrees of connected vertices as 
k n n(k) = ^2 k , k'P(k'\k), where P(k'\k) is the conditional 
probability that a given vertex k is connected to a vertex 



of degree k' . For a topology with no correlations among 
nodes' connectivity, the conditional probability as well 
as k nn {k) are independent of k. On the contrary, explicit 
dependence on k is a signature of non-trivial correlations 
among the nodes' connectivity and the possible presence 
of a hierarchical structure in the network topology. In 
the presence of correlations, a network can be classi- 
fied as one having "assortative mixing" or "disassorta- 
tive mixing", based on whether k nn {k) is an increasing 
or decreasing function of fc, respectively. In the former 
class, high degree nodes have tendency to be connected 
to nodes with large degrees, whereas in the latter class 
high degree nodes have a majority of low degree nodes 
as their neighbors. For a weighted network the weighted 
average nearest neighbors degree [7[ could be defined as, 

1 N 

3=1 

which measures the effective affinity to connect with 
high- or low-degree neighbors according to the magnitude 
of the actual interactions. The trend of fc™„ indicates the 
weighted assortative or disassortative nature of the net- 
work considering the strength of interactions among the 
nodes. 

We find (Fig. [5]) that ANI has uncorrelated structure 
for small degrees (k < 8) with no definite assortative or 
disassortative mixing. But for larger degrees (k > 8) 
ANI shows a clear disassortative mixing. This property 
of ANI is drastically different from that of WAN which 
was found to be having assortative mixing for small de- 
grees (k < 10) beyond which it approaches a constant 
value. The consistent disassortativity in ANI for higher 
degrees could be attributed to the political compulsions 
on regional and national hubs to provide connectivity to a 
large number of low-degree destinations. A hub in WAN 
is connected to many other hubs (in other countries) , thus 
compensating for the large number of connections with 
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low-degree airports. Since the number of hubs is limited 
in a national airport network, such as ANI, it could lead 
to disassortative topology. This feature of ANI is con- 
sistent with the nature of Airport Network of China @ 
which is reported to be disassortative. This could well be 
a feature of national infrastructure networks which have 
similar geo-political compulsions. The normalized corre- 
lation function (r) [H| (— 1 < r < 1) is another global 
quantitative measure of degree correlations in a network, 
which is zero for no correlations among nodes' connectiv- 
ity and positive or negative for assortative or disassorta- 
tive mixing, respectively. For ANI we find r — —0.4016, 
indicating disassortative mixing which is consistent with 
the observation in Fig. [6l Many technological networks 
and biological networks are found to be disassortative, 
while social networks are found to be assortative [23 |. 
Indian Railway Network Q has also been found to be 
disassortative. 

We find that weighted degree correlation, k™ n (k), is 
restricted in its range across k and is consistently higher 
than its unweighted counterpart, k nn {k), indicating to- 
ward bias in the traffic in ANI. The larger the degree 
of a node, the more pronounced is the difference. This 
implies that despite disassortative topology, the traffic is 
concentrated between high degree nodes. This is under- 
standable as most airline service providers concentrate 
on certain profitable, so called, trunk routes, thereby cre- 
ating high traffic corridors Q. 

VI. CONCLUSIONS 

We find that ANI, despite being small in size, has com- 
plex dynamics similar to those of bigger air transporta- 



tion networks. It represents an evolving transport infras- 
tructure of a developing nation. It presents a case of a na- 
tional airport network which has a different mechanism of 
formation than the global network. ANI, whose topology 
has a signature of hierarchy, has small- world network fea- 
tures and is characterized by a truncated scale-free degree 
distribution. Analysis of weighted ANI reveals clearer 
picture of the network dynamics. The traffic in ANI is 
found to be accumulated on interconnected groups of air- 
ports. It has disassortative mixing in contrast to WAN, 
with hubs having large number of low-degree neighbors. 
The topological disassortativity is offset by the traffic dy- 
namics as the traffic is concentrated between high-degree 
nodes. 

ANI is expected to grow at a rapid speed with ad- 
dition of airports, several low-cost air services Q, and 
importantly by increase in strength and complexity of 
interactions. It will be interesting to study the evolu- 
tion of this air transportation network, as it offers an 
excellent example to understand and model the topol- 
ogy, traffic dynamics and geo-political constraints which 
shape a national infrastructure of economic importance. 
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